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Abstract 

Despite several attempts, no efficient cluster algorithm has been con- 
structed for CP{N — 1) models in the standard Wilson formulation of lattice 
field theory. In fact, there is a no-go theorem that prevents the construction 
of an efficient Wolff-type embedding algorithm. In this paper, we construct an 
efficient cluster algorithm for ferromagnetic 5C/(A^)-symmetric quantum spin 
systems. Such systems provide a regularization for CP{N — 1) models in the 
framework of D-theory. We present detailed studies of the autocorrelations 
and find a dynamical critical exponent that is consistent with 2; = 0. 
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1 Introduction 



Cluster algorithms are among the most efficient tools for numerical simulations of 
nonperturbative lattice field theories and statistical mechanics models. The ffist 
cluster algorithm for discrete classical spins was developed by Swendsen and Wang 
for Ising and Potts models pOj. This algorithm was generalized by Wolff to 0{N) 
models formulated in terms of continuous classical spins [2]. In the 2-d 0(3) model 
the Wolff cluster algorithm has been extended to a meron-cluster algorithm which 
has been used to solve the complex action problem at non-zero 6'-vacuum angle 
CP{N — 1) models in two dimensions [H El El provide a generalization of the 
0(3) = CP{1) model alternative to 0{N) models. They have a global SU{N) 
symmetry and share several features with QCD, including asymptotic freedom, the 
nonperturbative dynamical generation of a mass-gap, a topological charge, and thus 
nontrivial 6'-vacua. Finding efficient algorithms for CP{N — 1) models is difficult 
in the standard Wilson formulation. In particular, a Wolff-type embedding cluster 
algorithm turned out to be inefficient for CP{N — 1) models [7J. In fact, there is 
even a no-go theorem that prevents the success of this type of algorithms IB]. Still, 
a rather efficient multigrid algorithm does exist (101 for ^ = 0. However, due to a 
very severe complex action problem, no efficient algorithm has ever been found for 
non-zero vacuum angle. Consequently, in the Wilson formulation numerical studies 
of 6*- vacua [II1II21 have been limited to moderate volumes, or are based on additional 
assumptions [Ullil 113113 • 

D-theory is an alternative formulation of field theory in which continuous classical 
fields arise dynamically as collective excitations by the dimensional reduction of 
discrete variables such as quantum spins |T71 UHl CHI- The discrete nature of the 
quantum variables allows the development of new types of cluster algorithms. The 
ffist cluster algorithm for quantum spins was developed in This algorithm 

works efficiently for 1-d ferro- and antiferromagnetic quantum spin chains, while the 
loop-cluster algorithm |2I1 123 works efficiently in any dimension. For example, the 
loop-algorithm has led to very accurate simulations of the 2-d Heisenberg quantum 
antiferromagnet 1221; which can even be performed in continuous Euclidean time 
The D-theory formulation of field theory has been used, for example, to construct 
an efficient cluster algorithm that solves the complex action problem in the 2-d 
0(3) model at non-zero chemical potential 123. In this paper we construct an 
efficient algorithm for the D-theory formulation of CP{N — 1) models. The D- 
theory approach allows efficient simulations both at ^ = and vr. In D-theory 
CP{N — 1) models are represented by ferro- or antiferromagnetic SU{N) quantum 
spin systems. Just like the SU{2) Heisenberg antiferromagnet, SU{N) quantum 
antiferromagnets can be simulated with the loop-cluster algorithm 123 which was 
recently reviewed in [2Z1I2H1- Here we construct an SU{N) loop-cluster algorithm 
for a ferromagnet to study CP{N — 1) models at 6' = 0, and we investigate their 
numerical efficiency. We find a dynamical exponent for critical slowing down which 
is consistent with z = 0. The D-theory method even solves the 6'-vacuum sign 
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problem, and has recently allowed us to confirm the conjectured first order phase 
transition [23 ED] at = vr [31] . 

The paper is organized as follows. We present the standard Wilson formulation 
of lattice CP{N — 1) models in section 2 and the D-theory formulation in section 
3. Section 4 describes the path integral formulation for CP{N — 1) models in terms 
of ferromagnetic quantum spins transforming in the fundamental representation of 
SU{N). The corresponding loop-cluster algorithm is constructed in section 5, and 
its efficiency is investigated in section 6. Finally, section 7 contains our conclusions. 



2 Standard Formulation of CP{N - 1) Models 

The manifold CP{N-1) = SU{N)/U{N-1) is a (2A^-2)-dimensional coset space 
relevant in the context of the spontaneous breakdown of an SU (N) symmetry to a 
U{N — 1) subgroup. In particular, in more than two space-time dimensions {d > 2) 
the corresponding Goldstone bosons are described by an x matrix-valued field 
P{x) e CP{N - 1) which obeys 

P{xf = P{x), P{xy = P{x), TrP(x) = 1. (2.1) 

For d = 2 the Hohenberg-Mermin-Wagner-Coleman theorem implies that the SU{N) 
symmetry cannot break spontaneously. Correspondingly, similar to 4-d non-Abelian 
gauge theories, the field P{x) develops a mass-gap nonperturbatively. Motivated by 
these observations, D'Adda, Di Vecchia, and Liischer |3] introduced CP{N — 1) 
models as interesting toy models for QCD. The corresponding Euclidean action is 
given by 

s[P] = I d^x ^Tiid^Pd^p] - leqypi (2.2) 

where g"^ is the dimensionless coupling constant. The topological charge 

Q[P] = ^Jd'x e,,TT[Pd,Pd,P] e U^iCPiN - 1)] = Z, (2.3) 

takes integer values in the second homotopy group of CP{N — 1)^, and 9 is the 
corresponding vacuum angle. Note that both the action and the topological charge 
are invariant under global transformations Q G SU{N) 

P{x)' = QP{x)Ql (2.4) 

The topological term explicitly breaks the charge conjugation symmetry C which 
acts as ^P{x) = P{x)*, such that Q[^P] = -Q[P]. It should be noted that C is 
not explicitly broken at = or vr. In particular, for 6 = n the Boltzmann weight 

"'^Note that the factor i in ea. H2.3|l is necessary for obtaining a real- valued quantity. 
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exp{i9Q[P]) = (— is C-invariant. For > 3 it has been conjectured that 
there is a first order phase transition at = tt, at which charge conjugation is 
broken spontaneously [21] ■ Recently, this has been confirmed numerically using a 
cluster algorithm similar to the one described below [SI]. For CP{1) = 0(3), on the 
other hand, charge conjugation remains unbroken, but there is still a second order 
phase transition at 6^ = vr |33 IS3 IS] • At that point the model corresponds to the 
k = 1 Wess-Zumino-Novikov-Witten conformal field theory [SHISSliSnj- 

It is straightforward to define CP{N — 1) models beyond perturbation theory in 
the framework of Wilson's lattice regularization. Then the field Px € CP{N — 1) is 
defined on the sites x of a quadratic lattice and the standard lattice action is given 
by 

^[^] = -;lETr[P.P.+A], (2.5) 
where fi is the unit-vector in the /i-direction. 



3 D-Theory Formulation of CP{N - 1) Models 

In this section we describe an alternative formulation of field theory in which the 
(i-dimensional CP{N — 1) model emerges from the dimensional reduction of discrete 
variables — in this case SU{N) quantum spins in d+1 space-time dimensions. The 
dimensional reduction of discrete variables is the key ingredient of D-theory, which 
provides an alternative nonperturbative lattice regularization of field theory. In D- 
theory we start from a ferromagnetic system of SU (N) quantum spins located at the 
sites x of a (i-dimensional periodic hypercubic lattice of size L'^. The SU (N) spins 
are represented by Hermitean operators = that generate the group SU{N) 
(e.g. the Gell-Mann matrices for the triplet representation of SU{3)) and thus they 
obey 

[T;, T^] = i5xyUcT:, Tr(T;T^^) = U^y^,,. (3.1) 

In principle, these generators can be taken in any irreducible representation of 
SU{N). However, not all representations lead to spontaneous symmetry breaking 
from SU{N) to U{N — 1) and thus to CP{N — 1) models. The Hamilton operator 
for an SU{N) ferromagnet takes the form 

H = -JJ2t:T^+v (3.2) 

x,i 

where J > is the exchange coupling. By construction, the Hamilton operator is 
invariant under the global SU{N) symmetry, i.e. it commutes with the total spin 
given by 

T^ = J:T^. (3.3) 
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The Hamiltonian H describes the evolution of the quantum spin system in an extra 
dimension of finite extent (3. In D-theory this extra dimension is not the Euchdean 
time of the target theory, which is part of the d-dimensional lattice. Instead, it is 
an additional compactified dimension which ultimately disappears via dimensional 
reduction. The quantum partition function 

Z = TTexp{-(3H) (3.4) 

(with the trace extending over the Hilbert space) gives rise to periodic boundary 
conditions in the extra dimension. 

For d > 2 the ground state of the quantum spin system has a broken global 
SU{N) symmetry. The choice of the SU{N) representation determines the sym- 
metry breaking pattern. We choose a totally symmetric SU (N) representation cor- 
responding to a Young tableau with a single row containing n boxes. It is easy 
to construct the ground states of the SU{N) ferromagnet, and one finds spon- 
taneous symmetry breaking from SU{N) to U{N — 1). Consequently, there are 
{N"^ — 1) — [N — 1)^ = 2N — 2 massless Goldstone bosons described by a field P{x) 
in the coset space SU{N)/U{N — 1) = CP{N — 1). In the leading order of chiral 
perturbation theory the Euclidean action for the Goldstone boson field is given by 

t 2ti 
S\P] = dt d'^x TilpsdaPd.P -— dr PdtPd^P]. (3.5) 
Jo J a'' Jo 

Here ps = Jn'^/A is the spin stiffness which is analogous to the pion decay constant 
in QCD. The second term in eq. ()3.5|) is a Wess-Zumino term which involves an 
integral over an interpolation parameter r. The point r = 1 corresponds to the 
physical space-time, such that P{x, t, 1) = P{x, t). At r = 0, on the other hand, the 
field takes the constant value P{x,t,0) = diag(l, 0, 0). At intermediate values 
of T the field is smoothly interpolated between the two limiting cases. There is 
no topological obstruction against such an interpolation because Ili[CP{N — 1)] is 
trivial. The integrand in the Wess-Zumino term is a total derivative. Hence, the 
integral depends only on the boundary values at r = 1, i.e. on the values of the field 
in the physical space-time. However, since H2[CP(A^ — 1)] = Z, there is an integer 
ambiguity of the integral depending on which particular interpolation is chosen for 
P{x, t, t). In order to make sure that this ambiguity does not affect the physics, the 
prefactor of the Wess-Zumino term must be quantized. Remarkably, this prefactor 
is just determined by the number of boxes n of the chosen SU (N) representation. 

Let us first consider the d = 2 case. For (3 = oo the system then has a sponta- 
neously broken global symmetry and thus massless Goldstone bosons. However, as 
soon as f3 becomes finite, due to the Hohenberg-Mermin-Wagner-Coleman theorem, 
the symmetry can no longer be broken, and, consequently, the Goldstone bosons 
pick up a small mass m nonperturbatively. As a result, the corresponding corre- 
lation length ^ = 1/m becomes finite and the SU{N) symmetry is restored over 
that length scale. The question arises if ^ is bigger or smaller than the extent /? of 
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the extra dimension. When ^ ^ (3 the Goldstone boson field is essentially constant 
along the extra dimension and the system undergoes dimensional reduction. Since 
the Wess-Zumino term contains dtP, it simply vanishes after dimensional reduction. 
Correspondingly, the action dimensionally reduces to 



which is just the action of the 2-d target CP{N — 1) model at 6' = 0. The coupling 
constant of the 2-d model is determined by the extent of the extra dimension and 
is given by 



Due to asymptotic freedom of the 2-d CP{N — 1) model, for small g"^ the correlation 
length is exponentially large, i.e. 



Here N/4tt is the 1-loop coefficient of the perturbative /5-function. Indeed, one sees 
that ^ 3> /3 as long as P itself is sufficiently large. In particular, somewhat counter- 
intuitively, dimensional reduction happens in the large f3 limit because ^ then grows 
exponentially. In D-theory one approaches the continuum limit not by varying a 
bare coupling constant but by increasing the extent f3 of the extra dimension. This 
mechanism of dimensional reduction of discrete variables is generic and occurs in 
all asymptotically free D-theory models [TTl UHl UHl ■ It should be noted that (just 
like in the standard approach) no fine-tuning is needed to approach the continuum 
limit. 

In order to show that D-theory indeed defines the same field theory in the contin- 
uum limit as the standard Wilson approach, we have evaluated a universal physical 
quantity in both regularizations jSl^. A convenient quantity is the universal finite- 
size scaling function F{y) = ^(2L)/^(L). Here C,{L) is the correlation length (ob- 
tained with the second moment method) in a finite system of size L, and y = ^{L) / L 
is a finite size scaling variable that measures the size of the system in physical units. 
Figure 1 shows Monte Carlo data for F{y) obtained both from D-theory and from 
the standard Wilson approach. Up to small scaling violations, the agreement of the 
two data sets confirms that, after dimensional reduction, the (2-|-l)-d SU{N) quan- 
tum ferromagnet indeed provides the correct continuum limit of the 2-d CP{N — 1) 
model. Thanks to the cluster algorithm, the D-theory framework allows calculations 
that are much more accurate than the ones using Wilson's approach. 

Let us also discuss the d > 2 case. Then dimensional reduction still occurs, but 
the Goldstone bosons remain massless even at finite f3. However, if f3 is decreased 
sufficiently, the SU{N) symmetry is restored and the Goldstone modes are replaced 
by massive excitations. By fine-tuning f3 to the critical value (again, as in the 
standard approach) one can reach a continuum limit of the target theory, provided 
that the corresponding phase transition is second order. 




(3.6) 




(3.7) 



^ oc exp{4:nPpjN). 



(3.8) 
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Figure 1: Monte Carlo data for the universal finite-size scaling function F{y) of the 
2-d CP{2) model. The filled circles represent D-theory data from a (2 + l)-d SU{3) 
quantum ferromagnet at [3 J = 6, while the open squares correspond to data obtained 
with the standard Wilson 2-d lattice field theory at = 2.25. 



4 Path Integral Representation of SU{N) Quan- 
tum Spin Systems 

Let us construct a path integral representation for the partition function Z of the 
SU{N) quantum spin ferromagnet introduced above. In an intermediate step we 
introduce a lattice in the Euclidean time direction, using a Trotter decomposition 
of the Hamiltonian. However, since we are dealing with discrete variables, the path 
integral is completely well-defined even in continuous Euclidean time. Also the 
cluster algorithm to be described in the following section can operate directly in the 
Euclidean time continuum [21]. Hence, the final results are completely independent 
of the Trotter decomposition. In d spatial dimensions (with an even extent L) we 
decompose the Hamilton operator into 2d terms 

H = Hi + H2 + ... + H2d, (4.1) 
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with 

Xieven Xiodd 

The individual contributions 

= -J T^T^^., (4.3) 

to a given if j commute with each other, but two different do not commute. Using 
the Trotter formula, the partition function then takes the form 

Z = lim Tr {cxp(-ei/i) cxp(-eii2)- exp{-eH2d)}^ . (4.4) 

M^oo 

We have introduced M Euclidean time-slices with e = j3/M being the lattice spacing 
in the Euclidean time direction. 

Inserting complete sets of spin states q e {u, d, s, ...} the partition function takes 
the form 

Z = ^exp(-%]). (4.5) 

[i] 

The sum extends over configurations [q] of spins q{x, t) on a (d + l)-dimensional 
space-time lattice of points {x,t). The Boltzmann factor 

M-l d 

exp(-5[g])=n n n exp{-s[q{x,t),q{x + i,t),q{x,t + l),q{x + i,t+l)]} 

p=0 i=l x=(xi,X2,---,x^) 

XiCven, t=2dp+i—l 

X Yl exp{-s[g(a;,t),g(a;-Fz,t),g(a;,t-M),g(a;-Fz,t-M)]},(4.6) 

x=(xi,X2,---,xj) 

Xiodd, t=2dp+d+i—l 

is a product of space-time plaquette contributions with 
exp(— u, u, u\) — exp(— s[d, d, d, d\) — 1, 
exp{—s[u,d,u,d]) — exp{—s[d,u,d,u\) — -[1 -|- exp(— eJ)], 

exp(— d, d, u]) — exp(— s[ci, u, u, d]) = -^[1 — exp(— eJ)]. (4.7) 

In these expressions the flavors u and d can be permuted to other values. All 
the other Boltzmann factors are zero, which implies several constraints on allowed 
configurations. 

In the following we will consider the uniform magnetization 

xi 

which represents one component of = J2x'^xj ^ ^^^^ ^ corresponding sus- 
ceptibility 

1 



where is the spatial volume of the system. 
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5 Cluster Algorithm for SU{N) Quantum Ferro- 
magnets 



Let us now discuss the cluster algorithm for the SU {N) quantum ferromagnet. Just 
like the original SU{2) loop-cluster algorithm jSHESj, the SU{N) cluster algorithm 
builds a closed loop connecting neighboring lattice points with the spin in the same 
quantum state, and then changes the state of all spins in that cluster to a different 
randomly chosen common value. To begin cluster growth, an initial lattice point 
{x, t) is picked at random. The spin located at that point participates in two pla- 
quette interactions, one before and one after t. One picks one interaction arbitrarily 
and considers the states of the other spins on that plaquette. One of the corners 
of this interaction plaquette will be the next point on the loop. For configurations 
Ci = [u, d, u, d] or [d, u, d, u] the next point is the time-like neighbor of (x, t) on 
the plaquette, while for configurations C2 = [u, d, d, u] or [d, u, u, d] the next point 
is the diagonal neighbor. If the states are all the same, i.e. for C3 = [u, u, u, u] or 
[d, d, d, d], with probability 

p=l[l + exp(-eJ)] (5.1) 

the next point on the loop is again the time-like neighbor, and with probability 
(1 — p) it is the diagonal neighbor. The next point on the loop belongs to another 
interaction plaquette on which the same process is repeated. In this way the loop 
grows until it finally closes. The cluster rules are consistent with detailed balance, 
i.e. 

pia)wiQ C,) = piC,)wiC, a), (5.2) 

where p{Ci) = exp(— s[Cj]) is the Boltzmann weight of the plaquette configuration 
and w{Ci Cj) is the transition probability to go from Ci to Cj. 

With this single-cluster algorithm, the uniform susceptibility of eq. ()4.9|l has an 
improved estimator 

X=(|C|), (5.3) 

given in terms of the cluster size 

|C| = ^ y 1. (5.4) 
2dM ^ 

This shows explicitly that the clusters are physical objects whose size is directly 
related to a physical quantity. 
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Figure 2: History of the magnetization of an SU (3) ferromagnet on a 320^ lattice at 
(3 J = 6.75 as a function of the Monte Carlo time (measured in units of sweeps with 
the multi-cluster algorithm). 

6 Efficiency of the Algorithm in the Continuum 
Limit 



In order to quantify the efficiency of our numerical method, we have used a multi- 
cluster algorithm for a (2 + l)-dimensional SU{3) ferromagnet in discrete Euclidean 
time to investigate the autocorrelation times of the uniform magnetization which 
gives the cleanest signal. A Monte Carlo history of this observable is shown in 
figure 2. This already indicates that autocorrelations are very much suppressed. 



L/a 


20 


40 


80 


160 


320 


640 




4.5 
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5.55 


6.2 


6.75 


7.45 




8.78(1) 


16.76(1) 


32.26(3) 


64.6(1) 


123.4(2) 


253(1) 


T 


0.721(3) 


0.719(4) 


0.724(6) 


0.719(1) 


0.715(2) 





Table 1: Numerical data for the correlation length C,{L) and the autocorrelation time 
T for an SU (3) -symmetric ferromagnet corresponding to a CP{2) model in D-theory. 
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We have analyzed the autocorrelation function for the magnetization which is 
illustrated in figure 3. In order to estimate the dynamical exponent z of critical 



0.1 



0.01 



0.001 



G O 


L/a = 


20 


□ □ 


L/a = 


40 


« O 


L/a = 


80 




L/a = 


160 




L/a = 


320 




Figure 3: Autocorrelation junction of the magnetization of an SU{3) ferromagnet 
as a function of the number of sweeps with the multi-cluster algorithm. In order to 
increase visibility, some points have been slightly displaced. 

slowing down, we have considered several lattices with decreasing lattice spacing, 
keeping the ratio y = ^{L)/L of the correlation length ^{L) (obtained with the 
continuous time cluster algorithm using the second moment method) and the size 
L fixed. Our data are summarized in table 1. Since the autocorrelation time r does 
not show any dependence on ^(L), critical slowing down is practically eliminated. 
In particular, our data are consistent with z = 



7 Conclusions 



We have constructed an efficient cluster algorithm for CP{N — 1) models within the 
theoretical framework of D-theory. The continuous classical CP{N — 1) field then 
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emerges dynamically from the dimensional reduction of discrete SU{N) quantum 
spins. Our cluster algorithm is a generalization of the loop algorithm for SU (2) 
quantum spins. We have investigated the algorithm's efficiency by estimating its 
autocorrelation times. Our Monte Carlo data show that critical slowing down is 
practically absent, i.e. the dynamical critical exponent is consistent with z = 0. 
This is thus another example where the D-theory formulation allows us to construct 
efficient numerical methods which are not available using the standard Wilson lattice 
field theory. In particular, the use of the discrete D-theory variables evades the no-go 
theorem that prevents the construction of efficient Wolff-type embedding algorithms 
for CP{N — 1) models in the standard Wilson formulation. Interestingly, using spin 
ladders with both ferro- and antiferromagnetic couplings, CP{N — 1) models can 
even be simulated at vacuum angle 6 = n [31j with no additional numerical effort. 
In particular, the sign problem that afflicts the standard approach is completely 
eliminated. Using the cluster algorithm described in this paper one can perform 
accurate numerical simulations of the CP{N — 1) mass gap for which no analytic 
results are available. Also the approach to the large limit can now be investigated 
in detail. It is also interesting to ask if D-theory allows us to construct efficient 
cluster algorithms for other field theories such as, e.g., SU{N) ® SU{N) chiral 
models with iV > 3 to which the no-go theorem also applies. Of course, the most 
challenging application of D-theory would be to non-Abelian gauge theories and, in 
particular, to QCD [HI CHI UHl • 
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